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Abstract. 

Simulations of a bilayer in a liquid solvent in the full range of areas from the floppy 
to the extended state also show the hole/tunnel spontanous formation. Whether this 
change of topology is preceded by topological fluctuations of the Gaussian curvature 
K{x, y) is an open question. Probability distributions of Gaussian curvature are trans- 
formed into those of the ordered inverse radii; these new shapes show a gradual change 
along the isotherm. The model of Max-Planck group with 5 beads is used, amended to 
represent amphiphiles in the spirit of the "coarse-grained" or "Martini" force fleld. 
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Introduction. 

Simulations of bilayers have greatly extended the knowledge about and the under- 
standing of these prototypes of membranes. Whereas most work was concentrated on 
the simulation of tensionless bilayers, the wider range of areas was examined as well, 
beginning with a small interval in the pioneering work of this kind[l], up to the entire 
range of the existence of a stable bilayer[2-7]. The simplest setup is the liquid system 
at constant temperature T, volume V, and constant N, the number of amphiphilic 
molecules, and A^^, the number of solvent molecules. With change of shape the area A 
varies and the lateral tension F with it; the "bilayer isotherm" r(^) results[2-7]. Its 
unusual and interesting features are briefly recalled in Section I. One may inquire as to 
the effect of temperature, but with one exception [7] and another [8] this has not been 
examined. We show some results elsewhere. With the increase in the area A at constant 
temperature, the stretched bilayer eventually gives way but it does so in an interesting 
manner; it does not break up in droplets or micelles but recedes to keep its surface 
density and creates a hole — a tunnel — through which the disjoined portions of solvent 
may connect. Thus the hole is filled with the solvent. The process of such sponta- 
neous formation, without any foreign body attempting to pass through the membrane, 
has not been studied much by simulation[9] ; we examine its appearing on the bilayer 
isotherm, in relation to the change of Gaussian curvature K. Normally this quantity is 
of no interest as it is null and constant, but with the hole formation the genus changes 
(from the value p = —2 appropriate for system with periodic boundary conditions) and 
therefore K changes. Section II is devoted to the Gaussian curvature. An interesting 
question is whether the sudden creation of hole is preceded by topological fluctuations, 
fluctuations in K. 

Section I. The bilayer isotherm and the model used. 

Fig.l shows a bilayer isotherm, the lateral tension F, plotted against the specific 
area per surfactant molecule a = 2A/N. The main portion of the graph shows the 
extended bilayer, (EX), where F is positive and increasing. Compressed bilayer reaches 
the tensionless state (which we call the "tensionless state nr.l"); at this area ai F 
vanishes. On further compression of the area the bilayer undergoes a transition to the 
"fioppy" state (FL). This is the small area of the graph where F is negative. Such 
transitions and the existence of the fioppy state are known from experimental studies 
of vesicles[10] and a counterpart is found here. The surprising feature is the constancy 
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of the negative T. All this can be understood and qualitatively described by the very 
simplest theory of undulations as put forward by Helfrich[ll]. We have also shown 
that the scaling with size is predicted by that theory and confirmed the prediction by 
our simulations of 4 systems of 4 sizes from N = 450 to N = 4000 [4]. Undulations 
of bilayers have been much studied, for example in[2-12], and the corresponding power 
spectrum of fluctuations, often under the name of the structure factor, S{q), has been 
determined repeatedly. S is expected to diverge on vanishing qf, — > 0+, as l/q^ for 
capillary waves and asl/q'^ for tensionless bilayer (in the "tensionless state nr 2"). Now, 
upon the transition to the floppy state the asymptote = is replaced by an asymptote 
at a positive value of q, q = qas > and S floppy (q) diverges as l/{q^ — q'^g). The region 
q e (0, qas) is not accessible. It can be shown[5] that such shift of the qf-asymptote is also 
predicted by the Helfrich undulation theory. The latter is a heuristic adaptation of the 
capillary-wave theory but it has been put on a firm footing in a mathematical paper[12] 
which derives all of it from a microscopic bending hamiltonian, by a coarse-graining 
procedure. It also follows from that theoretical description that the transition itself is 
not a true phase transition but a crossover - despite the appearances. 

On further compression of the already floppy bilayer there comes a point when 
the bilayer disintegrates into what appears to be a structureless foam and sometimes a 
"quad-layer". 

The picture described and seen in the left-hand side of Fig.l is valid for sufficiently 
large systems. For small bilayers the EX portion of the isotherm continues downwards 
to large negative values of F and then suddenly starts to increase strongly - so that 
the curve F(a) has a spike - and there is no flat section of constant, or nearly constant 
F[4,5]. In our experience at least 1200 amphiphiles were necessary to produce the 
expacted floppy plateau. 

The bilayer with the specific area increased above ai is in an extended state, gently 
undulating, with F increasing less than linearly, with a; then becomes stretched and 
finally does not disintegrate but recedes and creates a hole. The hole is nearly cylindrical 
and is filled with the solvent. The first published picture[3] of such a hole already showed 
that the bilayer recedes keeping its surface area; this was superseded by a more thorough 
investigation [9]. The lateral tension drops to much lower value. The downward trend of 
dV /da in Fig.l suggests that perhaps the hole is formed when dF/da = but on closer 
scrutiny the plot does not fully support that; the derivative seems to be still positive at 
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met ast able states. 

The details of the shape of a bilayer isotherm will depend on the model used and 
on the details onf the intermolecular force field. In Fig.l we show two sets of points, 
supplemented by a few points (shown with circles and line) at constant pressure and 
(slightly) variable density. This is meant to imply that it does not matter very much that 
we take isochoric data, at constant (overall) density. The main points (Diamonds) were 
obtained with the additional repulsive force (ARF) mimicking the hydrophobic effect 
after Goetz et al.[l]; the smaller set of data (Stars) - without that "Deus-ex-machina" 
addition. 

Fig. 2 shows a similar isotherm but for molecules with longer 8-bead tails. Since 
the system is also bigger, the floppy flat portion of r(a) is better developped and truly 
flat. 

The model we used in the past was introduced by the Max Planck group[l], as an 
extension of the very first modelling of amphiphilic molecules - as dimers in a solvent[13]. 
The liquid system is made of spherical particles, which are either solvent molecules or 
beads connected in linear chains of unbreakable (chemical) bonds to form amphiphilic 
molecules. In such a molecule, one terminal bead represents the polar head and the 
remaining - the nonpolar hydrocarbon chain. Here we use 5 beads {i.e. 1 for the head 
plus 4 for the tail) and exceptionally 9 i.e. one for the head and 8 for the linear chain 
as tail. The intermolecular energy is a sum of pair interactions, all of which are cut and 
shifted 6-12 Lennard-Jones functions: 

u{r)) =4e[{a/ry'' -{a/rf] (1) 

In general for 3 species (solvent "s", head "a", and tail-bead "t") the symmetrical matrix 
of parameters requires 6 energies e^/j and another requires 6 collision diameters aa/s- 
Practically in most if not all simulation work the collision diameters were taken equal for 
all pairs. Often the energy parameters were also made all equal[l,4,14] as well. Then, 
as one author noted[14] the molecule was not amphiphilic any more. 

Also, an extra term was added to substitute for the absent "hydrophobic force". 
It is well-known that a hydrocarbon molecule "prefers" its own environement to water 
and vice versa; it will experience an effective repulsion, in spite of all dispersion forces 
being attractive. Such an additional repulsive force (ARF) between the solvent "s" and 
the tail beads "t", has been added[l] and assumed to vary as with n = 8 or 10. 
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To improve on that model, we note that the essential feature of an amphiphilic 
molecule is the permanent chemical bonding of two antagonistic groups - a polar (hy- 
drophilic) head and the hydrophobic (non-polar) tail, most often a hydrocarbon tail. 
Accordingly, we treat these two groups of particles differently. The e's are taken all 
equal within the non-polar group 

^ = = ^at = ^sU (2) 
and augmented for the "polar group": 

^88 = ^sa = ^aa = (3) 

With u) > 1 the polar beads "s" and "a" are endowed with a stronger mutual attraction. 

In this simple manner the model is improved and the molecule cannot be said " not 
to be amphiphilic any more" . In most of our work we have used the value uj = 2. . For all 
intermolecular forces all cr's are taken equal. The cutoff at distance r = 2.5cr is common 
to the polar group. Next one manipulates further the potentials after Toxvaerd[15]. All 
potentials must be shifted if cut, so as to ensure continuity of forces. The s-t and a-t 
forces are made repulsive at all distances by changing the cutoff to r = a (and therefore 
the shift to +e). The intramolecular bonds are confined to the small neighbourhood of 
the optimal bond distance by a potential well with infinite barriers[4,15]. 

Our model can be viewed as a simplified version of the "Martini force field" [16], 
the latter also named a "coarse grained" model[17]. 

The Molecular Dynamics simulations were done in the canonical ensemble i.e. at 
constant T, V, N with Nose-Hoover thermostat. 

Section II. The hole formation and the topological fluctuations. 

Thw spontaneous hole formation in the stretched bilayer is a change of topology. 
It will affect the Gassian curvature K. Whether this sudden change of the surface (and 
of its gcmis) is preceded by fluctuations, e.g. by fluctuations of K, is not known. It is 
accepted that approaching even a sharp first-order order transition, "the system and its 
molecules exercise themselves" [18] in preparation to the transition. We have examined 
the behaviour of K along the bilayer isotherm reported above. If the Fourier coefficients 
have been determined, no matter how, one can construct the surface h[x, y) and from 
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the definition of K to obtain K{x^ y). Total, or overall, is then 



Kt = 




' / dxdyK{x, y) 
Jo 



(4) 



If a surface is given by the smooth function h{x,y) in the Monge gauge then 



K{x,y) = {h:cx * hyy - hly)/D^ 



(5) 



where 



D = 



(6) 



in the usual notation of first and second derivatives of h. From the Fourier analysis of 
the simulation data i.e. positions of heads, h is given as a finite sum 



where is a vector {27r /L){nx,ny) and R — {x,y). The g-vectors are usually chosen 
to fill a square centered at (0, 0) and a half of it is used in summation, coresponding to 
the square area 0<x<L, 0<|/<L. The other curvature, the mean curvature H is 
also obtained from h{x,y) by standard expression, as always with D ^ 1. 

Both were obtained from h{x,y) at a set of arbitrary values {xj,yj), chosen as a 
square grid sufficiently fine for trapezoidal integration. The latter produced the overall 
total Kt and total H'^, their dispersions, and a large set of values of K at {xj,yj) to 
be used for the construction of histogram. These computations were done along the 
entire simulation run and the sets of data obtained were used to construct histograms 



The histograms for H'^ are shown in Figure 3. The mean curvature squared H'^ 
decays monotonously attaining a rate not too far from an exponential decay in the 
latest stage of the decay, as can be seen from Figure 3 in semilogarithmic scale. The 
four points on the isotherm are specified in full detail below; the two upper curves for 
L = 53 and L = 52 are indistinguishable; that for 46 lies slightly higher than that for 
38. The next Figure, Fig.4, shows a portion of the same histograms, near H'^ = 0. The 
sequence is the same along the entire histograms (the curves never cross) and the Figure 
also shows the steep fall very close to H'^ = 0. 

The histograms for K are more interesting. P{K) appears as a Dirac delta-function, 
5(0), in agreement with the expectation that for a two-dimensional surface with periodic 



h{x,y) = ao/2 + ^ 



(In COs{qnRn) + 



(7) 



P{K)dK and P{H^)d{H^). 
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boundary condition the Euler genus equals -2 and therefore K = 0. Fig. 5 shows the 
remarkable distributions of K, in semilogarithmic scale to reveal the broadening of the 
(5-functions. 

The nature of the peak is examined in Fig. 6 by using a magnified scale; remarkably 
the discontinuity at = is approached with a finite slope and with an apparent lack 
of symmetry about K = 0. These two features: discontinuity with finite slope and 
no symmetry, are seen again in Fig.7. It shows the plots of the inverse 1/P{K), very 
convenient for the estimation of the slopes and extrapolation to the point of disconti- 
nuity. It is not clear if the extrapolation from the left produces the same value as the 
extrapolation from the right, in view of different slopes and second derivatives. The 
mathematical representation is an open question; a gaussian distribution is certainly 
excluded as there is no inflexion in sight. An inverse power is not acceptable either, as 
the limited plot of Fig.7 already shows. 

K was computed often, but not for flat bilayers; rather for those intricate surfaces 
formed in microemulsions and/or general minimal sur f aces (cf. [19,20]). In these nonpla- 
nar intricate shapes the Euler index is changed and is a most important quantity. 
There are no known expression for its distribution, even for one of the simplest mini- 
mal surfaces [19]. The shapes of those histograms, numerically computed [19], were very 
similar to those seen in our Figure 5. 

With the discontinuity a,t K = i.e. with the part with K < differing from the 
part with K > the task becomes even more complex. But this will be unnecessary as 
will become clear. 

All histograms as numerical data were obtained for different states of the bilayer, 
i.e. for different areas along the isotherm. Four points were chosen: (1) almost at 
the tensionless state, (2) in the middle of the isotherm, (3) near the breaking point of 
stretched bilayer and (4) at the breaking point where the hole-tunnel is formed. The 
details are as follows: T = 1.35, N — 1800,4 beads, p = 0.89 for 1 points; and then 
(1) L = 38.25, a = 1.62525; (2)L = 46.0, a = 2.3511; (3)L = 52.0, a = 3.004; and 
(4)L = 53.0, a = 3.1211. The last point is the point of hole formation; as mentioned 
earlier, there is hysteresis. There are metastable states when the hole is not formed for 
a very long time, for example, starting from the system at L = 52., after preliminary 
equilibration at L = 53., running with L = 53., we flnd no hole until 0.8E6 timesteps. 
Then expansion to L = 54. of the system with newly formed hole, stabilized the tunnel 
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structure, which on compression back to L = 53. did not disappear. The earher part of 
the data for L = 53., i.e. before the hole was starting to form, were included as (4) and 
is codenamed as "53h". As is apparent from Fig.l and 2, the hole formation makes the 
lateral tension F to fall to very different and low values [9]. 

For examining the approach to state (4) along the bilayer isotherm, i.e. in the 
sequence 1 — > 2 — > 3 — > 4 the plot of inverses 1/P{K) appears the best; the peak 
broadens as seen from Fig.7 but there is no regular pattern. 

The lack of symmetry about = is a disturbing feature. We turn therefore to 
the origin of these curvatures. 

The primary quantities are the principal radii of the surface h{x^ y) at that point, 
Ri{x^ y) and R2{x, y). The principal curvatures are defined as ci = l/i?i and C2 = 1/-R2- 
The principal radii are obtained as roots of a quadratic equation, usually written in 
terms of the c's 

a*x^ + b*x + c == a* {x — ci){x — C2) = 0. (8) 
As it turns out, H and K are related to the roots as 

iy = (l/2)(ci + C2) (9) 

K = C1C2 (10) 

Therefore H is simply related to one coefficient in the quadratic equation and K to the 
other. Simple algebra produces the solution, the roots, in terms of H and K, are 

ci = H- - K (11) 

C2 = i? + Vh^-K (12) 

In these relations of differential geometry H and K are just abbreviations for the com- 
binations of the first and second derivatives of h{x, y). Also 

Ri = (l/ci) = C2/K (13) 

R2 = (I/C2) = ci/K (14) 

A dump of the -R's and c's followed by the construction of their histograms, revealed 
interesting features. First, one radius is always very large in absolute value as compared 
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to the other. This is due to near-planarity of these surfaces. If |i?2| is 0(1), | i2i | is 10000 
or more. However, the signs vary; though ci < C2 always, not necessarily | ci |<| C2 | - 
so the smaller root may be the big one in absolute value and conversely. 

The output roots were therefore ordered according to their absolute values: if 
I ci |<| C2 I then Cb = ci, Cg = C2, and conversely, if | ci |>| C2 | then Cb = C2, 
Cs = ci. The convention is to use the subscript "s" for the radius smaller of the two, 
I Rb \>\ Rg I with Rg = 1/cs, Rb — l/cb and therefore for the larger of the two c's. It 
turns out that which root will produce Rb and which Rg, is decided by the sign of H, 
the mean curvature. Simply comparing the squares: > cf if < and conversely. 
And the value is 

(Rg)-^ = Cg = {\H\ +Vh^ - K) xH/\H\. (15) 

The histograms must be now transformed accordingly. First, histograms of unsorted 
roots Ci and C2 are shown in Fig. 8, for one point (2) (L = 46). A peculiar symmetry 
appears in that the branch for negative K is the mirror image of the other histogram 
for positive K: 

P(ci) = P(-C2) (16) 

The test for symmetry is better shown just for the positive semiaxis and with the inverses 
1/P(ci) and 1/P(c2), as in Fig.9. 

Now the histograms for sorted Cg and Cb take a very different look. Fig. 10 shows 
P{cb) and -P(cs), again for clarity for one state of the bilayer, namely (2) with L = 46. 
The larger radius Rb corresponds to Cb and P{cb)dcb has a spike at K = 0; the smaller 
radius is i?s = 1/cg and P{cs) vanishes at = and has a maximum. There is 
perfect symmetry now between the negative branch and the positive branch of the same 
histogram. P{cg) starts linearly from P(0) = 0) with the slope (1/4), goes through a 
maximum and then falls quasi-exponentially; in the final stage P ~ const * exp(— a * x). 
But as a fitting equation has very limited significance as it is only valid at very large 
values of the argument. Fig. 11 shows the pair from Fig. 10 now with semilogarithmic 
scale. It is seen that P{cg) decays almost exponentially, as exp(— ax) but with a = a{x) 
drifting along the curve to reach its final value at large values of the argument. 

The transformations to Rg and Rb substantiate and explain in simple terms the 
curious symmetries seen in Figure 8 but above all remove the lack of symmetries in 
histograms of K. The essential point is that it was not enough to use the principal 
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curvatures ci and C2; it is only the curvatures sorted according to their absolute values 
that had symmetrical histograms. It is physically convincing that the irrelevant big 
principal radius and its curvature C5, reflecting the overall near-planarity of the 
surface, were removed. The other and small radius Rs and its curvature Cs are the 
relevant quantitities. 

Now we compare -P(cs) for all state points (l)-(4); at last a pattern of regular change 
appears. As the area per surfactant molecule increases, Fig. 12 shows the sharpening of 
the distribution, a very small insignificant shift of the maximum, and constant initial 
slope at zero. The decay portion of P also changes regularly i.e. the decay is faster as 
the distribution gets sharper with the increase in area. 

From the histogram one obtains the average < Cg > and its dispersion (c^) — (cg)^; 
these are shown for all 4 state-points in Fig. 13. Just as the histograms P{cs)dcs there 
is smooth variation except for the additional point at L=38; this is a state point with 
ARF. 

In a related result, the variance of bare K is shown in Fig. 14. Taken as (KK) — 
{K){K), it shows good regular tends, unlike all untransformed histograms. Because 
K = C1C2, will be c^c| = c^c^ and apparently the regular variation of is not 
spoiled much by the presence of l/R^^^. Note that < K > may be slightly ofi' its true 
value of nil; wc take D ^ 1 so that K[x, y) does not vanish identically as it does in the 
small gradient aproximation D = 1. 

As mentioned earlier, it would be an unnecessary and tedious exercise to try to 
find a mathematicla representation fot the listograms of K. Also, even for one of the 
simplest minimal surfaces, this is not known[19]. Now, for the symmetric P{cs)dcs, it 
appears that a function 

P{x) ~ xexpl-ax""] {x > 0) (17) 

(with a mirror image for a; < 0) is not too far away, esp. for portions away from the 
maximmum of P. However, this attempt was not fully successful, as different least- 
square values for n were obtained: n = 0.668, n = 1. , and n = 1.12, 1.25. Trying 
to impose n = 1 for all was clearly not correct. And the maximum was not perfectly 
reproduced. Therefore a satisfactory mathematical representation, even for P{cs), was 
not found. 

These results have some further consequences of minor nature. The linear combi- 
nation ci -|- C2 or ci — C2 will be often swamped by the root bigger in absolute value and 
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therefore might not be very useful. We note that the quantity H'^—K = (l/4)(ci— 02)^ is 
always positive, so that the quadratic equation always has real roots and never complex 
roots. 

As the histograms of Cs were shown above, those for the radii, P{Rs)dRs, are not 
shown because their plots do not bring anything new. 

Summary and Discussion. 

We have found a new transformation of the curvatures H and the mean and the 
Gaussian, to quantities with exactly symmetric (about zero) histograms and varying 
smoothly along the bilayer isotherm. 

The numerical results were shown for 4 selected state points along the bilayer 
isotherm, one near the tensionless state no.l, one in the middle of the isotherm, one 
almost at the breaking point of hole formation, and one (metastable state) exactly at 
the hole formation. The data - positions of heads on both monolayers - were collected 
from simulation runs at T, V, N ^ one run producing e.g. 7600 data sets. 

The untransformed histograms - probability distributions - P{H'^)d{H'^) and 
P{K)dK, were examined and their plots are shown in the Figures. Their features led to 
reexamination of the derivation of these curvatures, as it is given in standard differential 
geometry; a new transformation resulted. In it, the principal radii are sorted according 
to their magnitude irrespectively of sign. This ordering is crucial; it is not sufficient to 
introduce principal curvatures ci and C2. 

These new objects, Cs = 1/Rs and Cb = 1/Rb, have symmetrical probability distri- 
butions which also vary smoothly along the isotherm. 

The resulting pattern of variation along the isotherm, can be interpreted by recalling 
that along the isotherm r(a), with the increase of the specific area a the bilayer is 
gradually stretched. Therefore predictably stretching impairs strong fluctuations of the 
monolayers and of the bilayer as a whole. Thus the distributions shown in Fig. 12 have 
smaller probability of large argument for the bilayer more strongly stretched. Similarly 
the dispersions fall with a, as shown by Fig. 13 and Fig. 14. 

There is no premonition of the breaking of the bilayer surface by creation of a hole 
- it appears to be a catastrophe without warning. As pointed out in the introduction, it 
is quite a special phenomenon, when the bilayer spontaneously recedes, trying to keep 
its structure and surface density, and makes room for a tunnel, the latter immediately 



11 



filled with solvent molecules. 
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Figure Captions. 

A short Caption is put under each Figure. 

In all Figures the scales are arbitrary with one exception of normalized distribution in 
Fig.l2. 

FIG.l. Bilayer isotherm T= 1.35, 4-1-1 beads, w/ extra force (Diamonds) and without 
(Stars); Circles - P = const. 

FIG. 2. Isotherm T=1.35 with (Diamonds) and without (Stars) the extra-force AFR; 
Boxes mark runs with holes. 8-1-1 beads. 

FIG.3. P{H'^)d{H'^) for L=53(hole),52,46,38. 

FIG.4. P{H^)d{H^) for L=53(hole),52,46,38. 

FIG.5. Histograms of Gaussian K for L=38,46,52,53,53h. 

FIG. 6. Selected peaks of Gaussian K; L=46(-|-),52(Stars),53h(Circles) and estd. 
maximum -line. 

FIG.7. Inverse of P(K) near peaks: L=38(+), 46(Stars), 53(Circles), 53h(Filled 
Circles). No pattern, no symmetry. 

FIG. 8. Principal curvatures; L=46. 

FIG. 9. Symmetries for principal curvatures. Proof that P{ci) = P(— C2). 

FIG. 10. Sorted roots (principal curvatures); the radius i?b is bigger in absolute 
value; = l/Rt - Plus signs, Cg = l/Rg - Stars. 

FIG. 11. Sorted roots (principal curvatures); the radius Rb is bigger in absolute 
value; Cb = l/Rb - Plus signs, Cs — l/Rg - Stars. 
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FIG. 12. P{cs) normalized to (1/2), with Rg being the smaller in absolute value; 
L = 38 (Boxes), L = 46 (Plus signs), L = 53 (Stars). 

FIG.13. Curvature c= 1/Rsm; (c) (Diamonds), (c^) (Stars), ((c- (c))^) (Boxes). 

FIG.14. Dispersions for H(Stars),H(sm.grad.) (Circles), K(Boxes). 
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FIG.1 .Bilayer isotherm T=1 .35,4+1 beads, w/ extra force(Diamonds) 
and without(Stars); Circles - P=const. 
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FIG. 2. Isotherm T=1 .35 with (Diamonds) and without(Stars) extra-force; 
Boxes mark runs with holes. 8+1 beads. 




FIG.3.P(H**2)d(H**2) for L=53(hole),52,46,38 




FIG.4.P(H**2)d(H**2) for L=53(hole),52,46,38 




FIG.5. Histograms of Gaussian K for L=38,46,52,53,53li 
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FIG. 6. Selected peaks of Gaussian K; 
L=46(+),52(Stars),53h(Gircles and estd. maximum -line) 
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FIG. 7. Inverse of P(K) near peaks:L=38(+),46(Stars), 
53(Circles), 53h(Filled Circles). No pattern, no symmetry. 
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FIG. 8. Principal curvatures; L=46 
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FIG. 9. Symmetries for principal curvatures. 
Proof that P(c-|) = P(-C2) and conversely 



0.035 



0.005 h 








1 


1 






+ 










+^ 








++ 








+ + 

+ + 





















■10 





c^=1/Rs c^=1/Rb 



FIG. 10. Sorted roots (principal curvatures); the radius R^j 
is bigger in absolute value; 

c^=1/Rb - Plus signs, c^=1/Rs - Stars 




FIG.1 1 .Sorted roots (principal curvatures); tine radius R^, 
is bigger in absolute value; 
Cg=1/Rs - Plus signs, C|3=1/R|3 - Stars 




FIG.12.P(c^) normalized to (1/2) with Rg being tine smaller 
in the absolute value;L=38(Boxes),46(Plus signs), 53(stars) 




FIG. 13. Curvature c=1/Rsm; average <c>(Diamonds), 
<cc>(Stars), and dispersion <(c-<c>)**2>(Boxes) 
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FIG. 14. Dispersions for H(Stars),Hsm.grad. (Circles), K(boxes) 



